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Abstract 

An influence of tlie anomalous Wtb coupling on forward-backward asymmetry in top quark 
pair production at the Tevatron is investigated taking into account decays of the top quarks 
to 6 fermion final states containing one charged lepton. To this end the most general effective 
Lagrangian of the Wtb interaction containing terms of dimension up to five is implemented into 
carlomat, a general purpose Monte Carlo program, which allows to compute automatically 
all necessary cross sections in the presence of anomalous vector and tensor form factors. A 
sample of results which illustrate little effect of the left- and right-handed tensor form factors on 
the tt invariant mass dependent forward-backward asymmetry and the charge-signed rapidity 
distribution of the lepton originating from the W boson from top quark decay is shown. 
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1 Introduction. 



The top quark is the heaviest particle ever observed, with mass close to the energy scale 
of the electroweak symmetry breaking. Therefore the top quark physics is an ideal place 
to look for non-standard effects which may reveal themselves through departures of the top 
quark properties and interactions from those predicted by the standard model (SM). The 
observation of a forward-backward asymmetry (FBA) in the top quark pair production in 
high energy proton-antiproton collisions at Tevatron [T], |2] that exceeds the SM expectation 
is an indication that this conjecture may be true. The CDF and DO Collaborations measured 
the total asymmetry Aft at a parton-level: 

AttiCBF) = 0.158 ± 0.075, Af(DO) = 0.196 ± 0.06, 

which is higher, but not inconsistent with the SM result. The asymmetry is zero in the lowest 
order of SM. A small asymmetry of A^ = 0.06 ± 0.01 arises at one loop QCD in the result of 
interferences of double-gluon corrections that differ under charge conjugation [3]. The CDF 
Collaboration finds that the asymmetry is a rising function of the tt invariant mass m^f, with 

Attima > 450 GeV/c^) = 0.475 ± 0.114, 

which is more than three standard deviations above the SM prediction in this mu region 
jl]. The DO Collaboration measured also a corrected asymmetry based on the lepton from 
a top quark decay to be 0.152 ± 0.040 which should be compared with the next-to-leading- 
order Monte Carlo generator result of 0.021 ± 0.001 [1]. Dedicated analyses of higher order 
contributions to the FBA of top quarks in the high invariant mass range of niu > 450 GeV/c^ 
show that the inclusion of the higher order QCD |5j and electroweak [6j, |^ corrections increases 
the one loop QCD prediction to some extent, but a 3cr deviation between the measurement 
and the SM prediction in this range still remains. Several new physics ideas, which alter the 
SM top quark production mechanism, have been invoked in order to explain the discrepancy 

El 

At the Tevatron, the top quarks are produced dominantly in pairs through the quark- ant iquark 
annihilation process 

qq ^ ti. (1) 

Creation of a top quark pair through the gluon-gluon fusion process, gg — )■ ti, that domi- 
nates the top quark production at the LHC, has much smaller cross section at the Teavatron. 
Moreover, it does not contribute to the FBA, as its initial state is symmetric under charge 
conjugation. Single top production processes, as e.g. qb — q't, qq' — tb or qg — )■ q'tb, have 
much smaller cross sections at the Tevatron, therefore their possible contribution to the FBA 
is neglected in the present work. 

Each of the top quarks of reaction ([1]) decays into a b quark and a W boson before hadronization 
takes place, and the W bosons decay into a fermion-antifermion pair each. The top quark 
pair production at the Teavatron is identified by selecting events where one W decays to qq' 
and the other to Wi. The experimental signature is an isolated electron or muon with large 
transverse momentum, a missing transverse momentum from the undetected neutrino and four 
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or more jets. At the parton level, one should consider reactions of the form 



uu{dd) — 7- bqq' blui, 



(2) 



where the quark q in the final state may be, but need not be identical with the initial state 
u or d quark. Any specific channel of ([2]) receives contributions typically from a few hundred 
Feynman diagrams, already at the lowest order of SM. For example, in the unitary gauge, 
assuming vanishing light fermion masses, = ma = rris = rrtf. = = 0, and neglecting 
the Cabibbo-Kobayashi-Maskawa (CKM) mixing between quarks, there are 718 lowest order 
Feynman diagrams for each of the reactions 



Examples of the Feynman diagrams of reaction ([3]) are shown in Fig. [H They include only six 
'signal' diagrams of ti production, three of which are depicted in Figs. [I](a) and[T]|^b) and the 
other three are obtained by permutation of identical u quarks. All the remaining diagrams 
constitute the off resonance background for the top quark pair production process. Some of 
them, as the one shown in Fig. [U^c), may contain a single top quark propagator, but most of 
the diagrams do not contain the internal top quark line at all, as the one shown in Fig. [11(d). 
Let us note that the Wtb coupling that is indicated by a black blob enters twice both in the 
tt production signal diagrams and the diagram with one top quark propagator. Obviously, it 
is not present in the off resonance background diagrams without internal top quark lines. 



Figure 1: Examples of the lowest order Feynman diagrams of reaction ([3]). Black blobs indicate 
the Wtb coupling. 

The presence of an anomalous Wtb coupling influences the top quark pair production in two 
basic ways. First, it changes the total decay width of the top quark, which substantially 
alters the total cross sections of any of reactions ([2]). Secondly, it changes the differential 
distributions of the final state particles, in particular of the final state lepton, which may have 
some influence on the tt production event reconstruction. 

Therefore, in the present Letter, the anomalous Wtb coupling of the most general form, with 
operators up to dimension five, is included in the theoretical analysis in order to see to which 
extent its possible modifications may change lowest order SM predictions for the ti invariant 
mass dependent FBA in the top quark pair production at the Tevatron. The question of 
whether the anomalous Wtb coupling may affect the asymmetry based on the charge and 
rapidity of the muon originating from the W boson from top quark decay will also be addressed. 



uu — i- bud bjj, 
dd — 7- ud bfi~Ufj_. 



(3) 
(4) 




b 
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2 An anomalous Wtb coupling 



The effective Lagrangian of tlie Wtb interaction containing operators of dimension four and 
five considered in tlie present Letter lias tlie following form [9]: 



'Wtb 



+ 



^ V 

71 

V* 



W-bY {fiPL + f^Pn) t - -^d^W-ba'^-' [f^PL + f^Pn) t 



mw 
1 



K^^' (fi^L + II'Pr) b d^W+ia^'' (J^Pl + J^Pr) b 



(5) 



where g is the weak coupling constant, rriw is the mass of the W boson, Pl = |(1 — 75) 
and Pr = |(1 + 75) are the left- and right-handed chirality projectors, a^'^ = | [7'^,7'^], Vtb 
is the element of the CKM matrix with the superscript * denoting complex conjugate, //', 



//^, //^ and 



R 



1,2, are form factors which can be complex in general. There are also 



other dimension five terms possible in Lagrangian ([S]) for off shell W bosons, but they have 
been neglected as they vanish if the W^s decay into massless fermions, which is a very good 
approximation for fermions lighter than the b quark. Therefore, in this approximation, Eq. ([5]) 
represents the most general effective Lagrangian of the Wtb interaction containing terms of 
dimension up to five. 

The lowest order SM Lagrangian of the Wtb interaction is reproduced by setting 



R 



m 
'2 



12 



(6) 



in ([5]). If CP is conserved then the following relationships between the form factors hold 



fR fh-' fLj 

Jl ~ Jl ; Jl ~ Jl 



f'. 



R' 



The Feynman rules resulting from 



cR 



and 4 independent form factors are left in Lagrangian 
Lagrangian ([5]) are as follows 




w:,q 



■pfj. 



riftPL + fi'p. 



R 



(8) 



and 




W-q 



_3_y* 



r {ftPL + I^Pr 



-^^^'^ [f^PL + f2^PB. 



(9) 



where g is a four momentum of the W boson outgoing from the Wtb vertex. 
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Direct Tevatron limits, that have been obtained by investigating two form factors at a time 
and assuming the other two at their SM values, are the following |12j^ : 



< 0.28. (10) 

The direct LHC limits that have been discussed in [T3] are still weaker. If CP is conserved 
then the right-handed vector coupling and tensor couplings can be indirectly constrained from 
the CLEO data on 6 — 57 [15] and from other rare B decays [16]. However, there is still some 
room left within which the anomalous form factors can be varied, in particular the tensor ones. 

The anomalous Wth couplings ([8]) and ([9]) are implemented into carlomat, a general purpose 
program for Monte Carlo (MC) computation of lowest order cross sections |10| . A new version 
of the program obtained in this way allows to make predictions for the top quark production 
and decay through different possible partonic subprocesses while taking into account complete 
sets of the lowest order Feynman diagrams and full information on spin correlations between 
the top quark and its decay products. The new version of carlomat can also be applied for 
studying anomalous effects in the top quark production and decay at the LHC, or in e^e~ 
collisions at a linear collider |17| . |18| . 



n 



R 



< 1.01, 



< 0.23, 



3 Results 



In this section, a sample of results that illustrate the influence of the tensor form factors of 
anomalous Wth couplings ([8]) and ([9]) on the ti invariant mass dependent asymmetry in the 
top quark production and on the charge-signed muon rapidity distribution at the high energy 
pp collisions at the Tevatron is shown. The results have been obtained with the current version 
of carlomat. 

The physical input parameters that are used in the computation are the following: the gauge 
boson masses and widths 

mw = 80.419 GeV, Tw = 2.12 GeV, mz = 91.1882 GeV, Tz = 2.4952 GeV, (11) 

the heavy quark masses and the Higgs boson mass 

mt = 172.5 GeV, rrift = 4.4 GeV, = 115 GeV (12) 

and the coupling constants 

aw = 1/132.5049458, a.(m^) = 0.118. (13) 

The QCD couplings are parametrized by = -yivraj. The electroweak coupling constants 
are parametrized in terms of g = ^/Anaw and the complex electroweak mixing parameter 
sin^ 9w = 1 — I M|, with the complex masses of the W and Z bosons My = my — imyTv, 

^After this work had been submitted for pubUcation new one-dimensional direct constraints at 95% C.L. on 
the form factors were announced by the DO Collaboration [T3]: |Vt&//^| < 0.93, |Vt6/^| < 0.13, |Vtfc/^| < 
0.06. 
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V = W, Z. The complex gauge boson masses together with the complex masses of the Higgs 



boson and top quark = — iniH^H and = yml — imtTt, where ReM^ > 0, replace 
masses in the corresponding propagators, both in the s- and t-channel Feynman diagrams. 
This choice of parametrizations is referred to as the 'complex mass scheme' in carlomat. The 
Higgs boson width is fixed at the lowest order SM value = 4.9657 MeV and the width of 
the top quark is calculated to the lowest order with effective Lagrangian 02]) for any specific 
choice of the form factors. 

The tt invariant mass dependent forward-backward asymmetry is defined by 

. . ^ _ (^i^y > 0, rritu) - a{Ay > 0, rriu^,) 
a{Ay > 0, mtt,i) + a[Ay > 0, ma^i) 

with Ay = Vt — Vt being a difference of rapidities of the t and i quarks with their invariant 
mass rriu within i-th bin. Since Ay is independent of boosts along the beam axis, asymmetry 
(1141) can be regarded as measured in the ti centre of mass system. 



The cross section of top quark pair production in pp collisions at ^/s = 1.96 TeV of Eq. ([T] 
is calculated by folding CTEQ6L parton distribution functions [19] with the cross section of 
hard scattering subprocess of the top quark pair production of the form ([2]), including all the 
subprocesses with uu and dd in the initial state and a single charged lepton in the final state, 
as e.g. processes (|3]) and (jlj). The factorization scale is assumed to be equal to a square of 
the reduced centre of mass system energy, s = X1X2S, with Xi {X2) being a fraction of energy 
carried by the initial state quark (antiquark). The ti production events are identified with 
the following acceptance cuts on the transverse momenta px, pseudoprapidities rj, missing 



transverse energy and separation ARi^ = \ {f]i ^ Vk)'^ + i^i ~ Vkf' in the pseudorapidity- 
azimuthal angle ((/?) plane between the objects i and k: 

PTi > 50 GeV/c, pTj > 50 GeV/c, \rii\ < 2.0, \rij\ < 2.5, 



> 20 GeV, ARu,ij,jj > 0.4. (15) 

The subscripts / and j in (ITSll stand for lepton and jet, a direction of the latter being identified 
with the direction of the corresponding quark. Cuts (fT5l) are rather restrictive. It has been 
checked by a direct computation that only events with the invariant masses of bqq' and bWi 
subsystems each close to survive. This means that the off resonance background contribu- 
tions are heavily suppressed and asymmetry (HM is in practice dominated by the events of ti 
production and decay. 

For the sake of simplicity, it is assumed that Vtb and form factors fP, Z/' and f^, i = 1,2 
of Lagrangian ([5]) are real. As the global fit combined with the SM constraints gives |Vtfe| = 
0.9991521o;qoqo45 |20| . a value of Va is fixed at Vtb = 1- Moreover, the vector form factors are 
assumed at their SM values of ([6]), = = 1, //^ = /f^ = and only the tensor form 
factors are being varied. 

In Fig. m asymmetry (fT4|) is plotted as a function of rriu, in bins of 50 GeV/c^ below 600 
GeV/c^ and 100 GeV/c^ above that. The plots in panels on the left hand side have been 
obtained with the complete set of the lowest order Feynman diagrams of each contributing 
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subprocess and those in panels on the right hand side have been obtained with the ti signal 
Feynman diagrams only. The result that corresponds to the form factors satisfying lowest order 
SM relations ([6]) is depicted with grey boxes in each panel, with solid error bars showing one 
standard deviation of the MC integration in separate bins. Boxes bounded by dashed lines 
show the asymmetry in the presence of two CP-even combinations of tensor form factors: 
f2 = I2 = 0.5, = /l^ = in the upper raw panels and = = 0.5, = = 
in the lower raw panels. The error bars drawn with dashed lines show the corresponding one 
standard deviation of the MC integration in separate bins. There are some fluctuations visible 
in separate bins, but they do not exceed 2a. Also the total lowest order asymmetry computed 
with carlomat is consistent with zero within one standard deviation of the MC integration, 
as it is shown in Table [H The asymmetry plots for other combinations of the tensor form 
factors, including the CP-odd ones, look very similar, so they are not shown here. Needless 
to say, the effect of the anomalous form factors becomes smaller if they are chosen within the 
recent DO limits |13| . 





Form factors 


mtt < 450 GeV/c2 


mtt > 450 GeV/c2 


Total 




= = f2^ = fl = 


0.09 ± 1.11 


-0.27 ±0.45 


-0.13 ±0.52 


f^ = 


fi = 0.5, fi = = 


0.07 ±1.23 


0.17 ±0.50 


0.13 ±0.54 


f^ = 


= 0.5, = fi = o 


0.39 ± 1.45 


-0.42 ±0.58 


-0.09 ±0.68 



Table 1: Asymmetry ( IT^ in % integrated in different ti invariant mass ranges. The vector 
form factors are fixed at their SM values of Eq. (jH]). The complete set of the lowest order 
Feynman diagrams is included. 

The charge-signed lepton rapidity distributions 1/cr da/d{qiyi), where yi is the rapidity of the 
charged lepton and qi is a sign of its electric charge, are plotted in Fig. [3] for different CP-even 
(upper raw) and CP-odd (upper raw) combinations of the tensor form factors of ([5]). The 
vector form factors are fixed at their SM values of Eq. ([6]). The complete set of the lowest 
order Feynman diagrams is included in the calculation. The SM result is shown in grey and 
the results obtained with non zero tensor form factors are depicted with boxes bounded by 
dashed lines in each panel. In spite of the fact that the anomalous Wtb coupling changes the 
total cross section of the top quark pair production by substantially altering the top quark 
width, the change in the charge-signed lepton rapidity distributions is hardly visible in the 
plots. 



4 Summary 

The tt invariant mass dependent FBA of top quark production and the charge-signed rapidity 
distribution of the lepton originating from the W boson from top quark decay at the Tevatron 
have been calculated to lowest order taking into account the anomalous Wtb coupling of the 
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Figure 2: Asymmetry (fT^ calculated with the complete set of the lowest order Feynman 
diagrams (left) and the tt signal diagrams only (right) for two different CP-even choices of 
the tensor form factors of (|5]). The vector form factors are fixed at their SM values of Eq. (Q. 
The asymmetry for all the form factors fixed by Eq. ([6]) is shown in grey in each panel. 



most general form, with operators up to dimension five [9]. It has been illustrated that even 
large values of the tensor form factors, exceeding the current limits [2], have rather little 
influence on the FBA. Also the charge-signed rapidity distribution of the lepton is very little 
affected by different CP-even and CP-odd combinations of the tensor form factors within the 
current limits. 
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